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Abstract 

Bekenstein's conjectured entropy bound for a system of linear size R and 
energy E, S < 2ttER, has counterexamples for many of the ways in which the 
"system," R, E, and S may be defined. Here new ways are proposed to define 
these quantities for arbitrary nongravitational quantum field theories in flat 
spacetime, such as defining R as the smallest radius outside of which only 
vacuum expectation values occur. Difficulties of extending these definitions 
to gravitational quantum and semiclassical theories are noted. 
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1 Introduction 



Bekenstein has conjectured that the entropy 5* of a system confined to radius 
R or less and energy E or less would obey the inequality (using units h = c = 

^3oltzmann 1) 

S < 2'kER. (1) 

He and colleagues have supported this conjecture with many arguments and exam- 
ples [1-18]. However, many counterarguments and counterexamples have also been 
noted [19-34]. Whether the conjectured bound ([1]) holds or not depends on what 
systems are considered and how i?, and S are defined. 

Perhaps the simplest procedure [H [10] would be to just consider quantum fields 
inside some bounded region within a sphere of radius R and put boundary conditions 
on the fields at the boundary of the region. However, this procedure leads to a large 
number of counterexamples to Bekenstein's conjectured bound. For example [20] , 
the Casimir effect can make E < Q for certain states of quantum fields confined 
within a certain regions of radius < i?, violating the bound. If states with E' < 
are excluded by definition, one can still consider a mixed state with arbitrarily small 
positive E that violates the bound. Even if E is redefined to be the nonnegative 
energy excess over that of the ground state [H [10], one can violate the bound by a 
mixed state that is almost entirely the ground state and a tiny incoherent mixture 
of excited states, at least if the entropy is defined to be S" = —trpliap [22]. If S is 
instead defined to be = In for a mixture of n orthogonal pure states (which would 
agree with 5* = —trpln p if the mixture had equal probabilities 1/n for each of those 
n pure states), then one can violate the bound by an equal mixture of the ground 
state and the first excited state of certain interacting fields with certain boundary 
conditions that have the two lowest states nearly degenerate in energy (separated by 
exponentially small tunneling effects) [20l [321 133]. If interacting fields are excluded 
from the definition of allowable systems, one can get a violation by considering a 
sufficiently large number N of identical free fields, giving n = N degenerate first 
excited states of finite energy but sufficiently large entropy S" = In n = In to violate 
([1]) [201 [21] • And even for a single free electromagnetic field, = Inn can exceed 
2'kER by an arbitrarily large factor by using boundary conditions corresponding to 
an arbitrarily large number of parallel perfectly conducting plates within the region 
of radius R [32], or by using boundary conditions corresponding to an arbitrarily 
long coaxial cable loop coiled up within the region [33] . 

However, other than in his papers with Schiffer [8], [10] , Bekenstein has generally 
advocating taking E to be the total energy of a complete system [H [3], [H [6], [151 [IB] • 
This would disallow using just the energy of fields within a bounded region with 
boundary conditions, since that would ignore the energy of whatever it is that is 
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providing the boundary conditions. Therefore, all of the counterexamples mentioned 
above would be excluded by this restriction. However, then the problem is to define 
what one means by the radius R of the system. In the weakly gravitating case 
(essentially quantum fields in flat Minkowski spacetime) that we shall focus on here, 
Bekenstein takes R to mean the radius of a sphere which circumscribes the system, 
which leaves the problem of what it means for a sphere to circumscribe the complete 
system. 

In quantum field theory in Minkowski spacetime, the complete system is the 
quantum state of the fields. Since the quantum fields extend all the way out to radial 
infinity, a sphere circumscribing the entire system would have to be at i? = cx), which 
makes the Bekenstein bound true (at least for states of positive energy and finite 
entropy) but trivial. To get a nontrivial bound, one needs to suppose that a sphere 
of finite R can circumscribe the system. For example, one might try to say that the 
sphere encloses all of the excitations of the fields from the vacuum. However, it is also 
hard to get this to occur for a finite R. For example, the wavefunction for any single 
particle state that is a superposition of energy eigenstates of bounded energy will 
not vanish outside any finite radius R, since a one-particle wavefunction that does 
vanish outside a finite region must be a superposition of arbitrarily large momentum 
components, which will have unbounded energy. Even if one looks at a composite 
system, such as a hydrogen atom, and ignores the fact that its center of mass will 
have amplitudes to be outside any finite sphere if it is made of purely bounded 
energy components, the wavefunction for the relative position of the electron and 
proton does not drop identically to zero outside any finite separation distance for 
states that are superpositions of energy eigenstates of bounded energy. In particular, 
even if one fixed the center of mass of a hydrogen atom in its ground state and 
ignored the infinite energy from the resulting infinite uncertainty of the center of 
mass momentum, the density matrix for the electron position would decay only 
exponentially with distance from the center of mass and never go to zero outside 
any sphere of finite radius R. 

Therefore, it is problematic to define the radius i? of a sphere circumscribing a 
complete system in any quantum field theory. This issue has not been addressed by 
Bekenstein and his collaborators, but without such a definition, there is no nontrivial 
formulation of the conjectured bound ([1]) for complete systems, but only its trivial 
truth for any complete system with positive energy and finite entropy that can only 
be circumscribed by the sphere enclosing all of space, R = oo. 

Here new ways are proposed to define systems and their radii R, energies E, and 
entropies S, so that for each, there is a bound on 5* for a given system as a function 
of finite R and E. These bounds will not have the form of Bekenstein's conjectured 
inequality ([I]), though in some cases they may obey that inequality. 
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2 Vacuum-Outside-R States 



The main new element of the present paper is a proposal is to define a system of 
radius R (in fiat spacetime for the present) not by imposing boundary conditions on 
the field itself, but by imposing conditions on the quantum state of the field so that 
outside a closed ball of radius R the quantum state is indistinguishable from the 
vacuum at some time. Such a state will be called a vacuum-outside- i? state. (For 
simplicity, set this time to be t = 0, and take the closed ball, say B, to be the region 
r < R on the t = hypersurface, where r is the standard radial polar coordinate 
giving the proper distance from the coordinate origin on that hypersurface.) In 
other words, a vacuum-outside- state of the system, say as expressed by its density 
matrix p, is such that the expectation value of any operator O which is completely 
confined to the region r > R when written in terms of field and conjugate operators 
at t = 0, is precisely the same as the expectation value of the same operator in the 
vacuum state |0><0|, 

tr{Op) =<0\O\0> . (2) 

In particular, all the n-point functions for the field and for its conjugate momentum 
in the state p are the same as in the vacuum state, if all of the n points are outside 
the ball of radius R and on the hypersurface t = 0. Of course, the ?7,-point functions 
need not be the same as their vacuum values if some or all of the points are inside 
the ball. 

If operators confined to the three-dimensional region r > R and t = (say C, to 
give a name to this achronal spacelike surface, the t = hypersurface with the central 
closed ball B, r < R, excluded) have the same expectation value in the vacuum- 
outside- i? state as in the vacuum state, the same will be true in any quantum field 
theory that I shall call "strongly causal" for all operators confined to the Cauchy 
development or domain of dependence [35] of C, the larger four-dimensional region 
r > R + \t\ (say D) that is the set of all points in the Minkowski spacetime such 
that every inextendible (endless) causal, or non-spacelike (everywhere timelike or 
lightlike), curve through such a point intersects the partial Cauchy surface C. Just 
as solutions of hyperbolic wave equations in D are determined by the data on C, 
so the part of the quantum state of a strongly causal field in D, as represented 
by the expectation values of operators confined to D, is determined by the part 
of the quantum state in C, as represented by the expectation values of operators 
confined to C. (For some interacting quantum field theories, the expectation values 
of operators confined to the three-dimensional spacelike surface C may be too ill- 
defined for these theories to be "strongly causal" in my sense, but a wider class of 
these theories may be "weakly causal" in the sense that sufficiently many operators 
smeared over, but confined to, an arbitrarily thin-in-time four-dimensional slab, say 
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E, containing C within D, have wcll-dcfincd expectation values that determine the 
expectation values of all operators smeared over, but confined to, any part of D.) 

Henceforth I shall restrict attention to strongly causal and weakly causal quan- 
tum field theories, calling them simply causal quantum field theories for short. I 
shall also assume, until discussing gravitational theories later, that any quantum 
field theory under consideration is a nongravitational Lorentz-invariant quantum 
field theory in Minkowski spacetime, and that it has a unique pure state of lowest 
Minkowski energy E = (the expectation value of the Hamiltonian H that gener- 
ates translations in the time coordinate t in some Lorentz frame, with the arbitrary 
constant in the Hamiltonian being adjusted to give the lowest energy state zero 
energy) . 

Therefore, for such a causal nongravitational quantum field theory in Minkowski 
spacetime, I shall propose that the radius R be defined so that all of the operators 
constructed from field and conjugate momentum operators smeared over regions 
confined to the region D, r > R + \t\ in some Lorentz frame, have in the particular 
quantum state being considered (a vacuum-outside- state) the same expectation 
values that they have in the vacuum state for that quantum field theory. The energy 
E of the state p can then be simply defined to be the expectation value. 



of the Hamiltonian H that generates time translations in the same Lorentz frame. 
Because the energy E has been defined to have the minimum value of zero for the 
unique pure vacuum state, there is no problem here with negative Casimir energies. 
In other words, the energy is that of the complete system over all of Minkowski 
spacetime. 

Obviously we would also like a definition of the entropy S that has a minimum 
value of zero, which it should attain for the pure vacuum state. One simple definition 
is the von Neumann entropy. 



using the density matrix p for the full state of the quantum field, over the entire 
Minkowski spacetime. 



3 Entropy Bounds for Vacuum- Out side-R States 



Now we may conjecture that for any vacuum-outside-/? state of any particular causal 
nongravitational quantum field theory in Minkowski spacetime, one which has the 
vacuum expectation values in the region D, r > R-\-\t\ (the region causally discon- 
nected from the ball r < R at t = 0), the von Neumann entropy is bounded above 



E = tr{Hp), 



(3) 




trp In p, 



(4) 
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by some function cTvn (depending on the quantum field theory in question) of the 
radius R and energy E: 

S,^<aMR,E). (5) 

Define this function a^^{R,E) to be the least upper bound on the von Neumann 
entropy of any state which is vacuum outside the radius R and which has energy E. 

In the case of a scale-invariant quantum field, such as a free massless field, or say 
a massless scalar field with a Xcj)'^ self-coupling potential, the least upper bound 
function a^-^{R,E) will actually be a function of the single dimensionless variable 

X = 27^RE, (6) 

say 

a^^{R,E) = cr^ix). (7) 

Bekenstein's conjectured entropy bound ([1]), if i?, E, and S were defined as done 
herein, would be a^]<f{R,E) < x, whether or not the quantum field theory is scale 
invariant, or 

5vN(i?, E) = = 2^^i^ < 1. (8) 

If the quantum field theory is scale invariant, we can define 

^ , , <yf<s(x) 

Bn X = 9 

X 

which should also be less than or equal to unity if Bekenstein's bound applies. 

For a set of one or more free massless fields and vacuum-outside- i? states with 
X ^ 1, one would expect that the highest entropy would be given by a mixed state 
that at t = is approximately a high-temperature {RT ^ 1) thermal radiation state 
for r < R, surrounded by vacuum for r > R. A high-temperature thermal radiation 
state has an energy density for massless fields of approximately a^T^, and hence an 
entropy density (4/3)arT^, where 

vr^, 7 , , , 

ar = —{nb + -nf) (10) 

is the radiation constant for ni, independent bosonic degrees of freedom for each 
momentum (e.g., Ub different spin or helicity states) and for nj fermionic degrees of 
freedom. Therefore, in this case with x ^ 1, 



X \ S^vr^x 



365x 



(rib + luf) 



111 



which is indeed less than 1, thus obeying Bekenstein's conjectured bound, for 

2V 2% 7 64nb + 56nf 

" ^ 3^ = 3^("^ + = 3645 ' ^''^ 
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if X is also large enough that Eq. (fTTj) is a good approximation. Thus one would 
expect that Bekenstein's conjectured bound, using the definitions above for _R, 
and S", holds for a fixed set of free massless quantum fields at sufficiently large 
X = 2nRE. 

On the other hand, the definitions above for R, E, and S still permit Bekenstein's 
conjectured bound applied to them to be violated for sufficiently small x, as we can 
see by the following construction: 

A way to construct vacuum-outside- states, quantum states of a free quantum 
field theory that have vacuum expectation values in the region D, r > R + \t\, 
is to apply to the vacuum state unitary operators constructed from fields and/or 
conjugate momenta smeared within the region r < R at t = 0. In particular, if h is 
an hermitian operator constructed from fields and/or conjugate momenta smeared 
within r < R at t = 0, then U = e*^ is such a unitary operator, and 

|^>= f/|0>= e*^|0> (13) 

is a pure quantum state that has precisely the vacuum expectation values in the 
region D. This result can be seen formally from the fact that any operator O 
confined to the region D (the four-dimensional region r > R + \t\) that is causally 
disconnected from the ball B (the three-dimensional region r < R on the t = 
hypersurface) commutes with the operators h and U that are confined to that 
hypersurface, [O, h] = [O, U] = 0, so 

tr{Op) =<iIj\0\iIj>=<0\U-^OU\0>=<0\U-^UO\0>=<0\0\0> . (14) 

If {hi] is a set of hermitian operators that each are confined to the ball B (i.e., 
are constructed from fields and momenta that are smeared only over that region), 
and if {g,} is a set of positive numbers that sum to unity, then 

p = ^g.e^'^'lOxOle-^'^' (15) 

i 

is a more general vacuum-outside- i? state, since this density matrix gives vacuum 
expectation values, tr{Op) =<0|O|0>, for any operator O confined to the region D 
that is causally disconnected from B (i.e., having no causal curves, either timelike 
or lightlike, intersecting both the ball B of r<R at t = and the region D with 
r>R+ \t\). 

Here let us consider the simple example in which i takes only the two values 1 
and 2, and hi = and /12 = h. Let qi = 1 — q and 52 = q, and let 

e*'^|0>= [/|0>= |?/^>= c|0> +s|l> (16) 
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in terms of a decomposition oi {ip > into the two orthonormal states |0 > and 
|1>= (|^> - <o|v^> \0>)/^l - |<0|V^>|2, so 



c=<0|?/^>=<0|[/|0>=<0|e^^|0>, (17) 



'1 - |<0|[/|0>|2 = ^1 - |c|2. (18) 

Note that |1 >< 1| by itself is not generically a vacuum-outside- i? state. 
Now Eq. f|T5l) gives the density matrix as 

p = (1 - g)|0><0| + glV^xV^I 

= (l-gs2)|0><0|+gcs|0><l| + gcs|l><0| + gs^|l><l|, (19) 

a density matrix in the two-dimensional space of pure states spanned by the two 
orthonormal pure states |0 > and |1 >. The two eigenvalues of this density matrix 
are, say, p and 1 — p (since their sum is trp = 1), with product 

y^pii-p) = \{[p+ii-p)r-[p'+ii-pm = l{Mp)r-Mp')]} 

= g(l-g)(l- |<0|[/|0>n = g(l-g)s^ (20) 

The expectation value of the energy of this mixed state is, since I have assumed 
H\0>= 0, 

E = tr{Hp) = q <^\H\^>= q <Q\U-^HU\Q> . (21) 

Then 

X = 27rRE = 27rRq <0\U'^HU\0> . (22) 
The von Neumann entropy of this mixed state is 

5'vn(i/) = —tr{p\np) = —p\np—{l—p)\n{l—p) 

= Ti ^ In — h Jl — Ay In -r- , 

i(i + v/r^) y V 1(1 + yr^) 

^ y[{l + y + 2y'')\n- + {l-\y-\y^)l (23) 

y 2 3 

a monotonically increasing function of y = q{l — g)s^ < 1/4, where the last approx- 
imate equality of Eq. (123|1 applies for ?/ <^ 1. 

As q and/or h is reduced toward zero, x, y, and S also decrease toward zero, but 
whereas x and y asymptotically decrease linearly with q, the dominant term of S 
has an extra logarithmic factor that grows with the reduction of y, so the ratio, 

B^^^^^yL- + l] (24) 
X 2-nRE x\ y ) ^ ' 

when ?/ -C 1, increases without limit as y is reduced toward zero. Therefore, when 
y is made sufficiently small (e.g., by making q sufficiently small), Bekenstein's con- 
jectured bound for the definition of i?, E, and 5* used here is violated. 
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4 Free Quantum Field Theory Examples 



Let us consider a specific example for the Hermitian operator h that is constructed 
from field operators confined to the ball B of radius R on an initial flat hyper- 
surface of Minkowski spacetime. Take the quantum field theory to be that of a 
single massless scalar field operator 4>. Consider the smeared linear Hermitian field 
operator 

X = J d^a:[F(x)0(t = O,x) + G(x)0(t = O,x)], (25) 

where F{x.) and G(x) are real functions of the spatial location x that are zero for 
|x| > R, so that X is made up of operators confined to the ball B, r < R at t — 0. 
Then for real parameters a and /3, let 

h^ax + Px", (26) 

which is thus also an Hermitian operator confined to the ball B. 

Then by expanding out 0(x) and H in terms of creation and annihilation opera- 
tors, one can show, after a certain amount of algebra that will not be repeated here, 
that 

c=<0|V^>=<0|t/|0>=<0|e*'^|0>= (1 - 2il3X)-^/^ exp ( (27) 

y 1 — 2ipX J 

^2 - I «iTTi„ lO - iio - /-. ..^o-r^O\_i/o / Ol ^ 



1 - |<0|[/|0>r ^ 1 - Ic^ = 1 - (1 + 4/?2x2)-V2exp (^^—^^j, (28) 
and 



<V^|i/|V^>=<0|[/-^ii't/|0>=<0|e-^'^i/e^''|0>= aV+4/3^XF = (a^X+4/?^X^)Z, 

(29) 

where 

X .<oix=|o>= /d3,,3^^w^W±Z5W:Z^ (30) 

J 47r''|x — yp 

and 

Y = XZ =<0\xHx\0>^ 2 / dH\F{^)\' + |VG(x)|^]. (31) 

Incidentally, I have not included individual higher powers of x h, because 
then expanding e*'^ into a power series in x and taking the expectation values gives 
divergent series when one uses the key intermediate results 

and 

^()\^rnrT^ni()^_ j mnim + n- 3)11 X^^'+'^y^Z m + neven , . 

1^ ^ I I m + nodd ■ ^ ^ 
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(The divergences arise from the rapid growth of the double factorials with their 
arguments. These double factorials arise from the counting of the number of pair- 
ings of the creation and annihilation operators in the powers of the x's and in the 
Hamiltonian H for the massless scalar field (p.) 

If one takes -F(x) and G'(x) to be spherically symmetric, say 



and 



F(x) = f 



G(x) = g 



|x| 

'r 

|x| 

'r 



R-'giu) 



(34) 



(35) 



with / and g being dimensionless functions of the dimensionless radial variable 
(hereafter to be called u or v) that vanish when the latter variable is greater than 
unity (corresponding to points outside the sphere of radius R), then 



X 



and 



du dv < 2uv In 



u + V 



u 



V 



f{u)f{v) + 



[u^ + v^) In 



u + V 



u — V 



2uv 



?7r /■! 

Y^XZ=- u'du[f{u) + g'\u)], 
R Jo 



g'{u)g'{v) 
(36) 

(37) 



where the prime on the function g denotes a derivative with respect to the argument 
(the dimensionless radius u = \:>c\/ R or v = |y|/i?). 

Now, if we take a density matrix of the form f|T9|) . let us try to maximize the 
product of the two nonzero eigenvalues of the density matrix. 



y = p{i-p) = g(l-g)s^ = g(l-g) 



-a'X 



1 + 4/52x2 



and hence maximize S^^{y) given by Eq. (!23l) . for fixed 

X = 2nRE = 2nRq <0\U'^HU\0>= 2nRZq{a'^X + A^^X^ 



(38) 



(39) 



Note that for fixed RZ = RY/X, the three quantities a, (3 (the coefficients of x and 
of in the hermitian operator h = ax + ^ind X =<0|x^|0> enter into this 

X and y only in the two nonnegative combinations a = o?X and h = 4/^2x2, and x 
depends only on g(a + 6). Then it is easy to see that for fixed q and fixed a + b, y 
decreases monotonically with b, so to maximize y and S'vn(|/) for fixed x, we should 
set P = in order to get b = 0, h = ax, 



a^X =<0\h^\0>. 



X 



2'KRZqa, 



(40) 
(41) 



10 



and 

y = g(l-g)(l-e-"). (42) 
Next, in our attempt to maximize y as a. function of x, we may continue to hold 

7 = 27rRZ (43) 

fixed and lience maximize y for fixed z = x/'y = aq. Then one can easily calculate 
that y = q(l — g)(l — e~^^'^) has its maximum at fixed z when 

a(l - -a), (44) 



2e«-2-a 2 ' 3 

giving 



and 



'jaie" - 1 - a) 1 2/-, 2 , 

x = 7ag = ^ -7a (1 a) (45) 

' ^ 2e'^ - 2 - a 2 ' ^ 3 ^ ^ ' 

y = \l-(l + a)e-"l f ~ ^ V ^ -0^(1 _ -a) ^ (46) 

where all the approximate equalities apply for a <^ 1. For fixed 7, a is given 
implicitly as a function of x by Eq. ( H5l) . and then inserting Eq. ( H6l) for ?/ into Eq. 
( !23|) for S'vN(y) gives the entropy (so far maximized over /3 and g) of the density 
matrix (JT9l) explicitly as a function of a =< 0|/;,^|0 > and hence implicitly as a 
function of x. In fact, for x -C 7, we get the asymptotic relation 

5,N(l/)~^(ln^ + l), (47) 

which of course exceeds x for sufficiently small x. 

As the final step in the maximization of the von Neumann entropy 5'vn(2/) of a 
density matrix of the particular form (fT9!) for fixed x = 2'kRE = 72; = 709, we note 
that maximizing y (and hence S^^{y)) for fixed x is equivalent to minimizing x for 
fixed y. Therefore, we need to minimize 7 = 27rRZ = 2ttRY/X. By looking at Eqs. 
( 130|) and ( l3Ti) for X and F, we see that the ratio Z = Y/X is invariant under any 
constant rescaling of the functions F{x.) and G(x) that appear as smearing functions 
for 0(t = 0, x) and (j){t = 0, x) in the defining Eq. ( |25l) for the linear hermitian field 
operators x and h = ax [now that we have set /3 = to drop the nonlinear term for 
h in Eq. fl26|) ]. The quantity 7 is also invariant under a rescaling of the radius R if 
F(x) and -F(x) depend only on x//2 and on some overall constant factor that can 
depend on R. 

Minimizing 7 = 2'kRY/ X is thus equivalent to maximizing X for fixed F, which 
by Eq. fl3T|) is half the integral of the sum of the squares of -F(x) and of the gradient 
of G(x). Because the double integral fl30l) for X is also quadratic in -F(x) and 
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in VG{x.) but has a positive-definite nonlocal kernel, maximizing X at fixed Y is 
best done with fairly smooth functions F{x.) and G{x). In particular, if -F(x) is 
expanded in spherical harmonics, one can readily see that the maximum is obtained 
by keeping only the spherically symmetric {£ = 0) terms. Also, since it is the dot 
product of V(j(x) and VG'(y) that enters into Eq. fl5Ul) . which generically dilutes 
its contribution relative to that of F(x)F(y) for the same values of the integrals of 
|F(x)p and of |VG'(x)p in Eq. ( 1311) . one can readily see that the maximum for X at 
fixed Y is obtained by setting G'(x) = 0, as well as choosing a spherically symmetric 
F(x) = /(|x|/i?)/i?2 as given by Eq. Then one gets that 
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J,'duJ,'dv2uv In ^J{u)f{v) 

One then sees that the minimum value for 7 is 

47r2 



(48) 
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A 



(49) 



where A is the largest eigenvalue of the weakly singular linear Fredholm integral 
equation of the third kind. 





U + V 


dv 2 In 


Jo 


u — V 



W[V) 



\w{u) 



(50) 



for < u < 1, with w{u) = uf{u) being the eigenfunction. 

For -F(x) = /(|x|/i?)/i?^ to be a smooth function of x, / should be a smooth even 
function of its argument u = |x|//2. This means that w{u) should be a smooth odd 
function of u, so we can expand it as an infinite sum of odd Legendre polynomials 

P2m-l{u), 



WiU) 



CmP2m-l{u). 



(51) 



m=l 



This expansion converts the integral eigenvalue Eq. fl50l) into the matrix eigenvalue 
equation 

00 00 



(52) 



n=l 



n=l 



where the matrix components are 
A 





U + V 


/ du c/f 2 In 


Jo Jo 


u — V 




2 



P2m-l{u)P2n-l{v) 



[1 — 4(m — n)^](m + n)(m + n 



and 



Br, 



du I dvP2m~liu)P2n-liv) 



Am — 1 



(53) 



(54) 
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[Actually, I cheated slightly in obtaining the explicit expression above for the 
matrix components Amn- I calculated An = 1 by hand, but when I tried to calcu- 
late the general A^n, I got finite sums that I did not readily see how to simplify. 
Therefore, I resorted to Maple. I did not quickly see how to get it to give me a 
simple general expression for Amn either, but in one afternoon I was able to get it 
to give me all the values for m < 10, n < 10 (45 different terms, since Amn = Anm)- 
The form of these terms was sufficiently simple that part way through their rather 
slow evaluation I was able to deduce the simple expression given in Eq. (l53l) . which 
indeed fit all 45 terms. So although I have not bothered to find a rigorous proof 
that Eq. ( l53l) is correct for all m and n not both smaller than 10, the fact that it is 
a very simple formula that works for all 45 smaller values strongly suggests that it 
is exact for all values of m and n. I could say that the proof is left as an exercise 
for the reader.] 

Maple readily solved the matrix eigenvalue Eq. fl52l) for various truncations of 
the infinite matrices Amn and Bmn- For example, 40-digit precision for 70 x 70, 
80 X 80, 90 X 90, 100 x 100, 110 x 110, and 200 x 200 truncations all gave the largest 
eigenvalue agreeing to 13 digits: 

A ^ 3.132010216749. (55) 

A 20-digit calculation of the 60 x 60 case gave the last digit 8 instead of 9 but was 
used to get the following approximate expansion of the eigenfunction corresponding 
to the largest eigenvalue: 

w{u) ^ + Pi{u) - 0.3968319408P3(m) + 0.0102661635P5(m) 

- 0.0070631 137P7(^i) - 0.0032552106P9(u) - 0.0018849293Pii(m) 

- 0.0011814233Pi3(m) - 0.0007878354Pi5(m) - 0.0005510316Pi7(m) 

- 0.0004002345Pi9(m) - 0.0002997366P2i(m) - 0.0002302203P23(m) 

- 0.0001806208P25(u) - 0.0001442924P27(m) - 0.0001170806Pi9(u) 

+ terms with coefficients less than 0.0001. (56) 

One can notice that only Pi{u) = u and Pziu) = —1.5m + 2.5m^ give large 
contributions to the eigenfunction, so one can get a fairly accurate estimate of the 
largest eigenfunction by taking even just the 2x2 truncation of the matrices, which 
gives the eigenvalue 



A2 = ^(^^ + ^) ^ 3.131921449343, (57) 
36 

which is smaller than the actual largest eigenvalue for the infinite matrices by less 
than one part in 35 283. An even simpler, but rather ad hoc, approximation is to 
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change the coefficient of the ^3(0;) term above to 0.4 and drop all the higher terms. 
Dividing this trial function for w{u) by u (and multiplying by 8 to avoid fractions 
in the answer) gives f{u) = 8 — 5m^, which may be inserted into Eqs. (1481) and fH^jl 
to give another estimate for A, 

1757 

Aest = -TT^r ~ 3.131907308378, (58) 
561 

which has almost 16% more error than A2, though this is still only a tiny error, 
being smaller than the actual largest eigenvalue for the infinite matrices by less 
than one part in 30434. Even the very crude constant trial function for f{u) gives 
an eigenvalue estimate, Acrudo = 3, that is smaller than the actual largest eigenvalue 
for the infinite matrices by only about 4.215%, or less than one part in 23. 

Using the approximation Maple gave for A, the largest eigenvalue of the infinite 
matrices, Eq. (H9ll then gives 

47r2 

7 = — ^ 12.604817632215, (59) 
A 

which can be used in Eq. f HTl) to get the asymptotic behavior of S^-^^x) at sufficiently 
small X. One can then see that this gives S^-s<f{x) > x for 

X < je^-^ ^ 0.000115, (60) 

or alternatively for 

5vN < 76^"^ ~ 0.000115. (61) 

Thus if the dimensionless energy, x = 27r RE, and the von Neumann entropy, S = 
S'vN = —trplnp, are sufficiently small, then with the definitions used here for these 
quantities, they can violate Bekenstein's conjectured entropy bound ([I]), 5' < x, 
though admittedly the range of x and 5* for which this happens is very narrow. 

We may now use the value of 7, given by Eq. (!59|) . in Eqs. (123|) . (145|) . and (H6|l 
to get B = S^-^/x as a precise implicit function purely of x, or, alternatively, to get 
both X and B as explicit functions of a =< 0\h'^\0 >. Of course, this is merely for 
one simple example of a one-parameter family of mixed states given by Eq. f|T9l) . 
with q given by Eq. (jHl) and {ip > given by Eq. ( |T6l) with (3 = so h = ax and with 
X given by Eq. ([25]) with G(x) = and Eq. §^ giving F(x) = f{\x\/R)/R^ with 
w{u) = uf{u) being an eigenvector corresponding to the largest eigenvalue. A, of the 
homogeneous linear Fredholm integral equation fl50l) . Therefore, it is not likely to 
give the maximum possible S as a function of x, which was called B^{x) in Eq. ([9]). 
However, this B{x) does give at least a lower bound on B^{x) for a single massless 
scalar field. 
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5 Conjectures for Entropy Bounds of Vacuum- 
Outside-R States 



I would conjecture that asymptotically at small x = 27iRE, the density matrix f|T9l) . 
with all the entropy-maximization procedures given above for a density matrix of this 
form, gives B{x) that does asymptotically approach the unknown global maximum 
function B-^{x) for a single massless scalar field. Therefore, if we divide x into the 
asymptotic form of S'vn for small x that is given by Eq. (117I) . I would conjecture 
that this gives the asymptotic form of the true upper bound, Bj<i{x), for very small 

X, 




with Eq. fISU]) giving 7 ^ 12.604817632215. One might expect a similar formula for 
other free massless fields, though perhaps each with a different value of 7. 

When X is not small, it is certainly not the case that the density matrix of fixed 
X needed to maximize the entropy is approximately of the simple rank-two form 
given by Eq. ([T9l) . One would surely need a more general vacuum-outside- state, 
with a density matrix obeying Eq. ([21) (giving vacuum expectation values for all 
operators not in causal contact with the ball r < R at t = 0), such as that given 
by Eq. (fT5l) . most likely with an infinite sum of terms and an infinite rank. I do 
not know how to proceed toward finding such a density matrix obeying Eq. ([2]) 
that would maximize B{x) = S/x at finite x that is neither asymptotically small 
or large. However, one might try using in Eq. ( |T5i) hi's that have the form given 
in Eq. (l26l) . with the F(x)'s and G'(x)'s of Eq. (p5|) being suitable eigenfunctions 
of the three-dimensional version of the integral equation ( |50ll . Even this wide class 
of examples may not be sufficient, since one could imagine instead constructing the 
Hermitian operators hi from smeared functions of the field 0(t = 0, x) and of its 
time-derivative (or conjugate momentum) 0(t = O,x) that are not merely linear as 
is the X given by Eq. (!25|) . Going to nonlinear Hermitian operators (other than the 
relatively simple considered above) leads to such a wealth of possibilities that I 
do not presently know how to proceed to obtain a true maximum for B{x) at fixed 
finite X, the postulated function B'^{x). 

In the more usual case of boundary conditions on the field, quantum states obey- 
ing these boundary conditions may be coherently superposed (i.e., the corresponding 
wavefunctions added, not just the density matrices) to get other states that also obey 
the boundary conditions. Then one can look for superpositions that diagonalize the 
Hamiltonian (i.e., energy eigenstates). From these, one can form a Gibbs ensemble 
to maximize the von Neumann entropy at a fixed expectation value of the energy. 

However, for the vacuum-outside-i? states considered here, coherent superposi- 
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tions of pure vacuum-outside- i? states are generically not vacuum-outside- i? states. 
(Of course, positive-weight combinations of vacuum-outside- i? density matrices are 
still vacuum-outside- i? density matrices when normalized, since the vacuum-outside- 
R condition, that all expectation values outside the ball r < R a.tt = are the same 
as the vacuum, is homogeneous and linear in the density matrix, though not in the 
wavef unction.) Therefore, the procedure for diagonalizing the Hamiltonian for such 
states fails. 

Indeed, one can see that none of the vacuum-outside- states, except for the 
vacuum itself, can be an energy eigenstate. This is because any state which is non- 
vacuum in a finite region at some time will inevitably have that region spread with 
time. For fields with linear field equations, the perturbations of the field itself will 
spread. But even for fields with self-coupling which allow classical field solitons 
that do not spread with time, any quantum state of the field which is non-vacuum 
at some time will inevitably have that region spread with time as a result of the 
quantum uncertainty principle. For example, suppose that there is some definition 
of the location and momentum of the soliton, such that the velocity of the location is 
proportional to the momentum. Then the position-momentum uncertainty principle 
will prevent one from having that the location remain, with certainty, within any 
finite region for an infinite amount of time; the quantum uncertainty of the position, 
if initially confined to a finite region, will inevitably spread to extend all over space. 
Therefore, the confined configuration cannot be stationary and hence cannot be an 
energy eigenstate. 

One can test my conjecture that Eq. fl62l) is the correct asymptotic form of the 
true upper bound on the entropy per x = 2ttRE by examining some other simple 
density matrices of the form given by Eq. (ITSl) that allow explicit evaluation of the 
von Neumann entropy. For example, the rank-three density matrix 

p = (1 - g)|0><0| + (g/2)e^°^|0><0|e~^°^ + (g/2)e-^"^|0><0|e^"'^ (63) 

has z = x/'j = aq, just like rank- two density matrix ( fT9l) when j3 = 0, and it has 
the three nonzero eigenvalues 

^ l-q{l-e-^) + l-q\l-e-^)^^l-qa+l-iq + q')a^ (64) 



P2 = lq{l-e-'n 



qa{l — a), (65) 
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and 

P3 = ^[1 - ^^(1 - e-'n] - - ^9(1 - e-2«)]2 - 2g(l - g)(l - e-)^ 

~ ^g(l-g)(l-e-'^)2^ig(l-g)a2, (66) 

where the approximate equations apply for very small a = o?X. If one chooses q to 
maximize the von Neumann entropy of this mixed state, 

'S'vN = -tr{p In p) = -pi Inpi - lnp2 - Vz Inps, (67) 

for very small a (equivalently, very small z), one gets that g ~ (1/3)(1 + 4a/3), 
so z = ag ^ iaj^i) + 4(a/3)^, which may be inverted to give a ~ 3z(l — 4z) and 
g ~ (1/3) (1 + 4^;). This then gives 




which has the same asymptotic form for small x as Eq. (1471) for the rank-two density 
matrix (1191) . 

Another example would be to consider the rank-five density matrix 

p=(l-g)|0><0| + (g/4) e*"^+*^>^' 1 >< 1 e-'°^-'^>^' 

+ (g/4)e*">^-*'^^'|0><0|e-^'"^+^'^^' 

+ (g/4)e-^°^+^^^'|0><0|e^"^-^^^' 

+ (g/4)e-'°^-'^^'|0><0|e*°^+*^^', (69) 

which gives z = xj'y = qia -|- 6), where, as above, a = o?X and h = 4/3^X^. Even 
though the eigenvalue equation is now a fifth-order polynomial equation, it appears 
that one may be able to use the symmetries of the problem to find the eigenvectors 
and eigenvalues explicitly, as functions of g, a, and 6, without requiring any roots 
higher than square roots. However, this seems to be messier than is worth doing 
here, so it shall be left as another exercise for the reader. 

Nevertheless, one can show that when a -C 1 and 6 -C 1, there is one eigenvalue 
near unity, one near ga, one near g&/2, and the remaining two are smaller by factors 
of the order of 2; = g(a + 6). Therefore, in this limit only the three largest eigenvalues 
contribute significantly to the von Neumann entropy, giving 

^vN ~ -(1 - ga - g6/2) In (1 - ga - g6/2) - gain (ga) - (g6/2) In (g6/2). (70) 

When this is maximized at fixed z = q{a+b), one finds that the first three eigenvalues 
need to be approximately in a geometric series (as, e.g., are all the eigenvalues of 
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the thermal density matrix for an harmonic oscillator), giving qh/2 ^ (ga)^(l + qa). 
Solving for qa and qb in terms of z and inserting this back into Eq. (!70l) gives the 
maximum von Neumann entropy 

^vN + J^-lnz] ~ - fln^ + l) (71) 

2 'J \ X J 

for the rank-five density matrix (169|) at fixed tiny z = x/7. This entropy is just 
slightly larger, by an amount roughly 2;^[ln(l/z) — (1/2)], than the corresponding 
maximum entropy (|68ll for the rank-three density matrix ( |63l) at tiny z = x/j, but 
it has the same asymptotic limit (given after the ~ sign). 

Therefore, although of course these three simple examples of finite-rank density 
matrices do not begin to exhaust the infinite set of possibilities, they give some 
support to the conjecture given above that Eq. ( !62l) is the correct asymptotic form of 
the upper bound for S/ {2nRE) when the denominator of this expression, x = 271 RE, 
is much smaller than unity. 

In the opposite limit, when x = 2ttRE is much greater than unity, we would 
expect that, at least for a scale- invariant field so that the energy E is large with 
respect to all relevant parameters with the same dimension (1/-R being the only 
relevant one if the field does not have a rest mass or other parameter setting a 
higher energy scale), the maximum entropy is given by Eq. ffTTl) for high-temperature 
thermal radiation, giving 

_^ 
x^ 



^n(x) ^ (72) 



with now 



■ 2^ 7 
365^ 8 ^' 



(73) 



no longer the (3 of Eq. ( l26l) that we have subsequently set to zero to maximize B{x) 
for the density matrix f|T9|) . It is tempting to combine this asymptotic formula for 
X ^ 1 with the asymptotic formula fl62|) conjectured above for x ^ 1 to conjecture 
that a reasonably good approximate formula for i?N( 

X^, clS db function of any x = 
2nRE, for a quantum field theory with a given set of massless fields, is 



5N(x)^-ln 1 + -^ , (74) 




where the constants (3 and 7 would depend upon the massless fields in the theory. 
For the single massless real scalar field that has been considered here, Uf, = 1 and 
n/ = 0, so Eq. (!26|) would give 




/J = ^ 0.364016115028, (75) 
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and Eq. (!59|) has already given 7 ^ 12.604817632215 for the single real massless 
scalar field. 

Another way to state this conjecture is to write 



B^{x) = -\n l + -^\Cix), (76) 




where C{x) is a correction factor yet to be found, and then conjecture that C{x) 
tends asymptotically to unity for both very small and very large x, and perhaps 
further to conjecture that C{x) is always relatively close to unity (e.g., say within a 
factor of two). This conjecture would then imply a conjectured entropy bound. 

An improved formula might be to write 

B^ix) = - In ^ ^ ^ ^^_i/4 jC'(x) (78) 

with B/C = (67)^/^ to fit the final 1 in the asymptotic formula ( !62l) for a; <^ 1, 
A-C = /37/4 to fit the asymptotic formula ([72]) for x > 1, and 2B-A^ + C^ = 7(5/2 
to fit the following two-term improvement to Eq. fl72|) : 



^N(a:)-^ + ^. (79) 

I have not tried to work out what 6 is. It would be straightforward to calculate, if 
it were the same as for the thermal state of a massless scalar field inside a sphere 
with Dirichlet boundary conditions on the field at the boundary r = R, but it is not 
obvious to me whether or not it is the same. 

Since Eq. f l78|) with the correction factor C{x) omitted should give a better 
asymptotic fit to B-^{x) than Eq. (CI]), I would expect that C{x) would generally 
be closer to unity than the corresponding correction factor C{x) of Eq. (176]) [not to 
be confused with the coefficient C in Eq. (178]) ]. But whether this is true over the 
entire infinite range of x remains to be seen. 

For free massive quantum fields, for fixed entropy one would expect that the 
energy would have to be higher, so an upper bound for a set of free massless quan- 
tum fields should also give an upper bound for a corresponding set of free massive 
quantum fields. Therefore, I would conjecture that for any given free quantum field 
theory, one can find a /3 and 7 (presumably with /3 obeying Eq. (175]) . and 7 some 
combination of eigenvalues of the appropriate integral equations) such that the in- 
equality (177]) holds with that value of /3 and 7. The conjecture might even be true for 
any reasonable interacting quantum field theory that is causal in the sense defined 
above, though then one might need a different value of f3. 
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6 Possibilities for Trying to Retain Bekenstein's 
Proposed Bound 



Returning to a consideration of how Bekenstein's proposed bound ([T]) fits with the 
results derived and conjectured here, one first notes that the results here violate 
([T]) for the von Neumann entropy of a vacuum-outside-i? mixed state with 
sufficiently small energy expectation value E, e.g., for x = 2'kRE < 0.000115 in the 
example above. However, even if one accepts the use of vacuum-outside-i? states 
for defining a finite size R, one might still object that the Bekenstein bound is not 
intended to be applied to the definition of E and/or 5" being used here. 

For example, Schiffer and Bekenstein [8] refer to "quantum states accessible to 
the field system with energy up to and including E." It could be objected that since 
the vacuum-outside- i? states considered above are not energy eigenstates, they are 
actually composed of states with energy both lower and higher than the energy 
expectation value that I have used as the definition of E. If one takes E to be the 
energy of one of the energy eigenstates that is sufficiently higher than the expectation 
value, then the Bekenstein bound ([T]) may be obeyed even in the examples I have 
given above that violate the bound when E is taken to be the energy expectation 
value. 

But if one takes this approach, it is hard to see how to give any content to 
the proposed bound for the vacuum-outside-i? states. Presumably not only is it 
the case that any vacuum-outside-i? state is not an energy eigenstate (since it is 
not stationary), but also it is surely the case that if any vacuum-outside- state is 
decomposed into energy eigenstates, it will include energy eigenvalues of arbitrarily 
large value. However, using an arbitrarily large value of E in the bound ([T]) makes it 
trivial, entropy less than or equal to infinity. Therefore, for the bound to have any 
content, we need to have a definition of E that gives finite values. The definition 
given by Eq. ([3]) above, E = tr{Hp), is surely the simplest, though others could be 
proposed. 

For example, one could propose instead that for a vacuum-outside- i? density 
matrix of the form (fT5|) . E could be defined as the maximum value of the expectation 
value of the Hamiltonian H in any of the normalized states e^'^'lO >< 0|e~*''' whose 
sum, weighted by the g^'s, forms p. However, using this definition would not avoid 
violations of Bekenstein's bound ([T]). For example, one could use the density matrix 
( !63|) with g = 1 so that the first term vanishes, and then the remaining two terms 
are of the form (fT5l) with hi = ax and /12 = — «X) and with qi = q2 = 1/2. Each of 
the two nonzero terms of the density matrix then gives the same energy expectation 
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value, and one can calculate that for small z = xj^ = a = o?X one gets 

~ z\\ - - (1 - z) Inz] ~ - fln^ + l) , (80) 
/ 7 V X / 

again violating Bekenstein's bound ([I]) for x < 0.000115. 

Yet another proposal would be that E be defined as -Emax? the maximum expec- 
tation value of the Hamiltonian in all of the eigenstates of the density matrix. In 
the example just discussed, 

p = p^\l><l\+P2\2><2\ (81) 
with pi = (1 + e~^°)/2, p2 = {1 — e"^'*)/2, and orthonormal eigenstates 

^ cos(ax)|0> (g2) 



and 

^^ sin(ax)|0> ^ (83) 
y/P2 



Then by using Eq. ( ]33ll . one can calculate not only that 

< 1 e"*'*! He'^' 1 > 1 e"^°^if e*"^ 1 > = < 1 e''^^ He'''^ 1 > 1 e^^^iJe"^"^ 1 > 

= a^Y = a^XZ = aZ (84) 

as given by Eq. fl29l) . but also 

< Ole-^'^iffe-^'^^ |0 >=< 0|e-*"^i/e-*"^|0 > = < Oje^^^iJe*^^ |0 >=< 0|e*"^iJe*°^|0 > 

= aZe-". (85) 

From these results and from the form of the orthonormal eigenstates given by Eqs. 
( I82l) and (J83l) . one readily obtains 

, , aZn - e-") 1 . 

1 + e / 

and 

ry 

H22=<2\H\2>= — (87) 

Then if one takes -Emax, the larger of Hn and if22, namely iJ22, as the definition of 
E, one sees that it has the positive lower limit Z = 7/(27ri?) as one takes a (and 
hence Svn) to zero, so with this definition one does not get a violation of Bekenstein's 
bound ([1]) in this example. In particular, for this example 

Be =— = i^^{-l(l + e-2-)ln[-(l + e-2-)] 

- ^(l-e-^'^)ln[^(l-e-^")]}<l. (88) 
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Whether this definition of Be^^^ always gives a result less than unity for all vacuum- 
outside-i? states, thus agreeing with Bekenstein's bound, remains to be proven, 
but the extremely meagre evidence that I have does not seem to contradict this 
conjecture. 

Another objection that might be made against the violations of Bekenstein's 
bound ([1]) using vacuum-outside- i? states to define i?, using the expectation value 
of the Hamiltonian given by Eq. ([3]) as the definition of the energy and using 
the von Neumann entropy given by Eq. (jlj) as the definition of the entropy is 
to demand that the entropy instead be given by a microcanonical ensemble rather 
than by Eq. (jll) applied to any mixed state. In other words, instead of allowing a 
generic vacuum-outside-i? mixed state or density matrix p obeying Eq. ([2]) for all 
operators O completely confined to the region D, r > R+ \t\, one might propose that 
Bekenstein's conjectured bound should only be applied to density matrices made up 
of equal mixtures of n orthogonal vacuum-outside- i? pure states. (These are rank-n 
density matrices with precisely n nonzero eigenvalues, all equal to and with the 
corresponding set of n orthonormal eigenvectors all being vacuum-outside-i? pure 
states.) The entropy of such a density matrix whose nontrivial part is proportional 
to the identity matrix in the n nontrivial dimensions is then 5* = Inn. 

Again, I do not have evidence that Bekenstein's conjectured bound ([T]) is violated 
for such restricted vacuum-outside- i? density matrices. However, it is a rather severe 
limitation to restrict the discussion to such a small subset of vacuum-outside-i? 
density matrices, a subset of measure zero in the space of all such density matrices. 
Furthermore, it appears rather difficult to find many explicit examples of precisely 
orthogonal vacuum-outside- i? pure states. 

For example, for each fixed choice of the two functions F{x.) and G(x) in Eq. (!25l) . 
the constants a and P in Eq. (126|1 give a two-parameter family of vacuum-outside- i? 
pure states of the form e*"^^*^-^^ |0 >, and then Eq. (1271) and its trivial generalization 
gives the inner product between any two states among this two-parameter family. 
However, none of these inner products are zero for finite a's and /3's (and hence for 
finite expectation values of the energy), so none of these states are orthogonal for 
fixed -F(x) and G{x.). 

Of course, if one combines various ones of these nonorthogonal pure state density 
matrices to get a mixed density matrix and then finds the eigenvectors of that 
density matrix, they will form an orthonormal set of density matrices, such as the 
set |1 >< 1| and |2 >< 2| of Eqs. (!82l) and (1831) . However, these density matrices 
are not by themselves vacuum-outside- states, but only when they are combined 
with the particular eigenvalues pi and p2 given just before Eq. fl82l) . Hence they 
cannot be used in a different linear combination (e.g., with pi = P2 = 1/2) to get 
a vacuum-outside-i? state that is an equal-weight combination of n orthonormal 
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pure-state vacuum-outside- i? states. 

The only explicit vacuum-outside- i? state orthogonal to the vacuum itself (the 
trivial vacuum-outside- i? state) that I have found so far is the pure state 

|V'>=exp(me"'^^')|0> (89) 

with real a. and /? (not the same a and /5 used elsewhere in this paper) chosen so 
that 

<0|^>=5]^^(l + 2n/3X)-^/2^0. (90) 

n=0 

Using Maple, I found a numerical solution at 

a ^ 4.727048274, p ^ 1.536994796/X (91) 

I have not worked out the expectation value of the energy, <ip\H\ip>, of this pure 
state, but I suspect that it is greater than ln2/(7r_R), so that the entropy of an equal 
mixture of this state and of the vacuum. In 2, would be less than 27ri? times the 
expectation value of the energy in this mixed state (half the expectation value of 
the energy of the pure state \ip >< ipl, since the vacuum half of the mixed state 
contributes zero to the expectation value of the energy). If so, then this example 
would not be a counterexample to Bekenstein's conjectured bound ([1]) restricted 
to microcanonical ensembles that are equal mixtures of orthogonal pure vacuum- 
outside- i? states. 

It would be interesting to find the lowest-energy vacuum-outside-i? state or- 
thogonal to the vacuum itself and see whether its energy is indeed not more than 
ln2/(7ri?), but I do not see how to do this at present. More generally, one would 
like to find, for each positive integer n, the set of n mutually orthogonal vacuum- 
outside- states (possibly, but not necessarily, including the vacuum itself) such 
that the sum of the n energy expectation values, say Es, is minimized. Then if one 
finds that Eg > (n In n) /(27ri?) for each n, then Bekenstein's conjectured bound ([1]) 
will be obeyed for these microcanonical ensembles of vacuum-outside- i? states. 

One way to look for other pure vacuum-outside- i? states orthogonal to the vac- 
uum would be to choose some non-hermitian operator, say n, that is confined to the 
ball r < i? at t = 0, and consider the one-complex-parameter (two-real-parameter) 
set of states 

\^{C;K)>=e'^^-'^''^\0> (92) 

for the complex parameter C. For a generic such k, < 0\iIj{C; k) > would be a 
complex function of C (not analytic, since both C and its complex conjugate C 
appear in the definition of \ip{C] k) >), and a simple parameter-counting argument 
suggests that there should be discrete complex values of C at which <0\ip{C; k) >= 
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0, giving a pure vacuum-outside- i? states orthogonal to the vacuum, though of course 
for particular k's, the number of such discrete solutions for C may be zero. One 
might extend this method to try to find n mutually orthogonal vacuum-outside- i? 
states; this would require n different operators /tj and n(n+l)/2 complex parameters. 
The obvious problem for carrying out this procedure explicitly is that for most sets 
of operators Kj, the inner products (functions of the complex parameters) would be 
difficult to evaluate. 

Perhaps a compromise to the stringent requirement of a microcanonical ensemble 
of n equally- weighted orthogonal pure vacuum-outside- i? density matrices is simply 
to use the von Neumann entropy S'vn = —trplnp, which equals Inn for a micro- 
canonical ensemble, but require that it be at least as large as In 2, the minimum 
nontrivial value for a microcanonical ensemble. Then one might conjecture that the 
bound ([1]) is correct for vacuum-outside- i? states such that S = S'vn > 1^2. 

Alternatively, one might replace Bekenstein's conjectured bound ([1]) with the 
weaker conjectured bound 

S <27rER + \n2, (93) 

still using S = S'vn = —trplnp, E = tr{Hp), and restricting to vacuum-outside- i? 
density matrices p obeying triOp) =<0|O|0> for all operators O totally confined 
to the region D, r > R+ \t\, that is not in causal contact with the ball B, r < R 
at t = (no causal curves connecting these two regions). ( Equivalent ly, one may 
require that triOp) =<0|O|0> for all operators O that commute with all operators 
defined totally on the ball B.) 

7 Other Ways to Define a Radius R 

So far I have been considering only the new proposal to define R by restricting to 
vacuum-outside- -R states. However, one might ask whether there are other ways to 
define a radius R for a class of states for which one is seeking a bound on the entropy 
S' as a function of R and of the energy E. 

One proposal that is very close to my proposal of vacuum-outside- i? states is a 
proposal for what might be called stressless-outside-i? states, states such that on the 
t = fiat hyperplane of the Minkowski spacetime that I have always been assuming 
so far in this paper, the expectation value of the regularized stress-energy tensor 
operator, T^^, is zero everywhere outside the radius R (as it is everywhere for the 
vacuum state), 

V(x) = tr{T^,{t = 0,^)p) = (94) 
for all |x| > R. Alternatively, one might restrict to what might be called energyless- 
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outside- i? states, 

e(x)=roo = tr(roo(t = 0,x)p) = (95) 

for all |x| > R, for which the expectation value e{x) of the regularized energy density 
operator, Tqq, at t = vanishes outside the radius |x| = R. 

Of course, all vacuum-outside- i? states are also stressless-outside-i? states, and all 
stressless-outside-i? states are also energyless-outside-i? states, but I do not know 
whether the converses of these statements are true. If they are not both true, 
there would exist energyless-outside-i? states, and possibly also stressless-outside-i? 
states, that are not also vacuum-outside- i? states. If there is indeed a broader class 
of states than vacuum-outside- i? states, whether stressless-outside-i? states and/or 
energyless-outside-i? states, then the corresponding entropy maximization function 
(TvnI-R, E) would be expected to be larger for the broader class of states. 

One can try to define R for even broader classes of states, not by requiring 
that some expectation values vanish for r > i? at t = 0, but instead by using the 
spatial distribution of some quantity to define an effective radius R. One way that 
first comes to mind is to use some spatially-dependent real weight function iy(x) 
coming from the quantum state to define R as an rms value of r = |x|: 

Of course, the weight function should be such that both the numerator and the 
denominator are finite and have the same sign (which without loss of generality will 
be assumed to be positive), at least for the class of states to be considered. 

An obvious simple choice of the weight function is the energy density expectation 
value e(x). Then the denominator of Eq. (!96l) for i?^ is the total energy, which 
is positive for a nontrivial state. However, the numerator is not positive for all 
nontrivial states, as one can see from the following argument: Motivated by the 
state with locally negative energy density given by Kuo and Ford [37], consider the 
state 

|^>= a|0> +/?|2>, (97) 
where |ap -|- = 1 and |2> is a two-quantum state of energy 



E2=<2\H\2>= J d-^x <2\Too\2> (98) 
and with mode functions that are sufficiently localized that 

RIE2 = f d^x <2|Too|2> (99) 



is finite. Because Tqq is the regularized (e.g., normal-ordered) energy density oper- 
ator, < 0|Too|0 >= and / d^x < 2|Too|0 >= / d^x < 0|Too|2 >= 0. For a generic 
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two-quantum state |2>, 

C = J d^x <2\Toq\0> (100) 
will be a nonzero complex number. Then Eq. fl96|) for Ty(x) = e{x.) gives 

Since can be an arbitrary complex number, R^, can take any real value if C 7^ 0, 
including zero and negative values. Even if one restricted to states for which R^r 
is positive, this quantity can be made arbitrarily small, and then if such a state 
has finite energy and positive entropy (e.g., by being a mixture of the vacuum state 
|0><0| and of {ipxtpl), it can make B = S/{27tRE) arbitrarily large. 

If one did want to use £:(x) as the weight function iy(x) in Eq. fl96|l . one would 
have to restrict the states so that i?^ cannot be too small for states of finite energy 
and nonzero entropy. One way that might work would be to restrict the states 
to those in which £:(x) is nonnegative everywhere, unlike the state \ip >< ipl for 
sufficiently large —Ca/p. 

Another option that might work for all sufficiently localized states would be to 
choose a weight function W{x.) that is nonnegative for all states. Examples of this 

would be EJ=oELo(v)'> i^P'^ i^u'rp', {rl^rl^rPy, {rl^rlr^rlf, 

(tr(0(t = O,x)p))2, (tr(0(t = O,x)p))2, (tr(: ^^(t = 0, x) : p))^, (tr(: ^^(t = 0, x) : p))^, 
(tr(: 0'^(t = 0, x)0;^(t = 0, x) : p))^, etc., and positive powers of these positive func- 
tions of x at t = 0. 

The fourth quantity above is the square of r/fr^, which itself usually seems to 
be positive everywhere, though it might be some restriction of states for this to be 
true everywhere for all states in the class. 

If one used one of these positive quantities, or one of the usually positive quanti- 
ties just for states in which it is everywhere positive, as a weight function Vr(x) for 
defining R by Eq. ( l96l) . one would again presumably get some upper bound on the 
entropy S as a function of R and of E (depending on how one defined S and E and 
what further restrictions one puts on the states). However, I have not investigated 
what these relations might be. 



8 Difficulties with Entropy Bounds in Quantum 
and Semiclassical Gravity 

So far I have been restricting attention to nongravitational quantum field theories in 
flat Minkowski spacetime. However, since the original motivation for Bekenstein's 
conjectured entropy bound came from quantum considerations of gravitating black 
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holes, it is interesting to consider whether a similar entropy bound can be applied 
in quantum gravity. 

Here I must admit that I see serious problems in attempting to apply the bound to 
quantum gravity. Assuming that the quantum part of quantum gravity is sufficiently 
similar to the ordinary quantum theory of nongravitational systems, the entropy S 
might still be a well-defined quantity, at least for a complete system, such as the 
entire universe (though if the ultimate quantum gravity theory specifies a unique 
quantum state, there may be no option as to what the entropy is) . The energy E is 
more problematic, at least if the universe is not asymptotically fiat, though if one 
can restrict to quantum states in which the universe is asymptotically fiat, then E 
also might have a good definition in quantum gravity. However, what I don't see 
how to give a good precise definition for is the size R. 

The main problem is that states of quantum gravity should be coordinate invari- 
ant, so it is hard to see how to say that some state is confined to a radius R or has 
this size. How would one define the center with respect to which the state is within 
a distance Rl Furthermore, if the state is an asymptotically flat one with energy E, 
the gravitational field of this energy should extend all the way out to spatial infinity, 
so in that sense it seems that the state cannot be confined to be within radius R or 
have vacuum properties outside that radius. 

It is not that I have a rigorous proof that an entropy bound such as Eq. ([T]) 
cannot be applied in quantum gravity, but I just don't see how it can be applied. 

The situation seems somewhat more hopeful in semiclassical gravity, in which one 
has quantum field theory for nongravitational fields on a classical curved spacetime 
(perhaps whose Einstein tensor is proportional to the regularized expectation value 
of the stress-energy tensor of the nongravitational quantum fields). In this case one 
can imagine defining the equivalent of vacuum-outside- i? states of energy E in the 
following way: 

Take an asymptotically flat spherically symmetric spacetime which has a totally 
geodesic Cauchy hypersurface (with zero extrinsic curvature, say at t = 0), about 
which it has time-reflection symmetry {t — > —t). Use a Schwarzschildean radial 
coordinate r, the circumference/ (27r) of each symmetrical sphere. Outside the r = R 
sphere on this t = hypersurface, assume that the spatial metric and the expectation 
value of all operators confined to this region are the same as that of the static 
spherically symmetric asymptotically Schwarzschild semiclassical metric with ADM 
mass E and quantum state that is the semiclassical version of the zero-temperature 
Boulware state for this metric. The entropy S can then be the von Neumann entropy 
= —trplnp of the quantum state of the nongravitational quantum field in the 
classical curved metric. 

If one applies this definition to all possible states of this form, then it seems 
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one can easily violate a bound of the form ([I]) by a state with arbitrarily large 
entropy by having the r = R sphere be a neck separating the asymptotically flat 
exterior with an interior on the t = hypersurface that is almost an entire three- 
sphere of arbitrarily large size filled with thermal radiation. In other words, take 
the interior of the r = R two-sphere to be the moment of maximum expansion of an 
almost-complete large k = 1 radiation Friedman- Robert son- Walker model, and take 
the exterior to be a moment of time-symmetry of a nearly empty approximately 
Schwarzschild metric. If the interior three-sphere radius is a ^ R, giving an interior 
volume going as a^, the semiclassical Einstein equations imply that the radiation 
energy density must go as in Planck units at the moment of maximum expansion, 
so the temperature T goes as a ' and the entropy density goes as or as a ' . 
When this is multiplied by the volume, one gets an entropy going as a^/^, which can 
be made arbitrarily large for fixed R and E (the asymptotic ADM mass) by making 
the interior size a arbitrarily large. 

One might seek to avoid this violation of ([T]) by restricting the states to which it 
is conjectured to apply to exclude this example of a huge interior universe separated 
from an asymptotically fiat exterior by a relatively small neck. 

One way to do that would be to demand that inside the r = R two-sphere on 
the hypersurface of time symmetry, there are no round two-spheres of radius greater 
than R (i.e., topological two-spheres with intrinsic two-metrics that are those of the 
standard unit round two-sphere multiplied by a constant that is larger than R^). 
This would exclude an interior that is approximately a large round three-sphere 
of radius a ^ R, since such an interior region would have round two-spheres of 
radii r ^ a. However, it still appears to allow a very long throat of radius near R, 
which by its arbitrarily great length could have arbitrarily large volume and hence 
arbitrarily large entropy S for fixed R and E. 

Another way to restrict the states so that they might possibly obey a bound 
similar to ([1]) is to demand that the evolution of the semiclassical geometry give 
a nonsingular metric over the whole of R'*. This would exclude the examples of a 
large internal approximate three-sphere and also the long internal throat, since these 
examples would be expected to collapse gravitationally to singularities. Only in cases 
in which the metric is not too much different from flat spacetime would one expect 
that no singularities develop from gravitational collapse, and in these cases one might 
expect an entropy bound not too different from its fiat spacetime form. However, for 
the restricted states to be sufficiently broad to encompass most of the semiclassical 
gravity generalizations of the allowed states (e.g., vacuum-outside- i? states) in the 
nongravitational theory, the semiclassical Einstein equations should give nonsingular 
evolution in these cases of sufficiently weak gravity. Since the semiclassical Einstein 
equations are of higher order in time than the ordinary Einstein equations with a 
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classical source, it is not clear that this will be the case, and so one might need to 
find a suitable semiclassical gravity theory before trying to apply entropy bounds. 

9 Conclusions and Acknowledgments 

In conclusion, we have found that one can formulate precise definitions for entropy 
bounds of a complete quantum field system (i.e., one not restricted to the interior 
of some boundary) by giving precise definitions for the size R of the system, at least 
when the metric is classical so that sizes can be unambiguously defined. In particu- 
lar, R may be defined for vacuum-outside- i? states as the largest round two-sphere 
on a suitable t = hypersurface, outside of which all of the operators have the same 
expectation values as in a suitable vacuum state (e.g., the ordinary vacuum state for 
nongravitational fields in Minkowski spacetime, or a Boulware-type quantum state 
in semiclassical gravity). Other values of R may also be defined, such as the rms 
value of r with a suitable weight function dependent upon the quantum state of 
the field. On the other hand, for a fully quantum gravity theory, it appears to be 
difficult to give an unambiguous definition of a size i? of a system, so it is not clear 
there how to define a bound for the entropy S in terms of the energy E and a size 
R. 

Discussions with Valeri Frolov, Jonathan Oppenheim, and L. Sriramkumar have 
been helpful. This work was supported in part by the Natural Sciences and Engi- 
neering Research Council of Canada. 
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